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A  model  is  presented  for  approximating  load-diffusion  from  axially  loaded  fibers  embedded 


in  elastic  matrices.  The  fundamental  elastostatic  solutions  used  are  for  a  point  force  and  a  point 


dilatation  in  either  a  fully-infinite  or  semi-infinite  space.  Tangential  tractions  across  the  fiber-matrix 
interface  are  included  explicitly  in  the  analysis.  The  model  is  applied  to  the  three-dimensional 
analogs  of  Melan's  first  problem  and  Reissner's  problem  and  comparisons  are  made  with  exact 


results  in  the  case  of  the  former  to  help  establish  the  validity  of  the  model. 


INTRODUCTION 


"v  The  ability  to  analyze  load-transfer  in  fiber-matrix  systems  which  are  illustrative  of  those 
that  exist  in  fiber-reinforced  materials  is  fundamental  to  the  study  of  how  such  materials  behave  in 
application.  Our  ability  at  present,  however,  to  rigorously  solve  such  problems  in  the  realm  of 
three-dimensional  elasticity  is  limited  to  a  few  isolated  results  involving  infinite  fibers  bounded 
along  their  entire  length  to  fully-infinite  matrices.  Two  noteworthy  examples  are  the  load-diffusion 
from  an  axially  loaded  fiber  [1]  and  the  load-absorption  by  a  broken  fiber  in  a  remotely  stressed 
medium  [10].  A  much  more  interesting  class  of  fiber-matrix  system  involves  fibers  embedded  in 
semi-infinite  half-spaces.  This  type  of  problem  is  pertinent  to  the  study  of  how  fiber-bridging  in 
the  wake  of  crack  advance  serves  to  isolate  the  crack  tip  from  applied  far-field  loadings  and  thus 
increases  the  fracture  toughness  of  a  material.  The  analytic  complexity  of  such  problems,  however, 
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discourages  any  attempt  at  a  rigorous  solution.  What  is  being  proposed  in  this  paper  is  a  model  for 
approximating  load-diffusion  in  these  systems. 

One  such  model  has  already  been  developed  by  Muki  and  Sternberg  [1,2,3]  and  used  to 
study  such  problems  as  load-transfer  to  a  half-space  from  a  partially  embedded  axially  loaded  rod 
[2]  and  load-absorption  by  a  semi-infinite  fiber  in  a  remotely  stressed,  fully-infinite  matrix  [3]. 
Muki  and  Sternberg's  model  replaces  the  fiber-matrix  system  of  the  problem  with  an  extended 
matrix  occupying  the  volume  originally  containing  both  the  fiber  and  the  matrix  and  possessing  the 
same  elastic  properties  as  the  original  matrix.  This  extended  matrix  is  in  tum  reinforced  by  a 
"fictitious  stiffener"  whose  modulus  of  elasticity  when  taken  in  sum  with  that  of  the  extended 
matrix  is  equal  to  that  of  the  original  fiber.  This  stiffener  is  taken  to  be  a  one-dimensional  elastic 
continuum  bonded  to  the  extended  matrix  in  such  a  way  that  the  axial  strain  in  the  stiffener  is  equal 
to  the  average  extensional  strain  of  the  extended  matrix  in  the  volume  occupied  by  and  in  the 
direction  of  the  original  fiber.  Poisson's  effect  in  the  stiffener,  and  therefore  in  the  fiber,  is  not 
taken  into  account.  Finally,  "bond-forces"  are  regarded  as  body  forces  uniformly  distributed  over 
disks  perpendicular  to  the  axis  of  the  fiber  and  the  load  carried  by  the  original  fiber  is  equated  with 
the  sum  of  the  stiffener  load  and  the  resultant  load  carried  by  the  extended  matrix  in  the  bonded 
region. 

A  variation  of  Muki  and  Sternberg's  model  was  used  by  Pak  in  a  study  of  flexure  of 
partially  embedded  fibers  under  lateral  loads  [4],  The  concept  of  a  "fictitious  stiffener"  replacing 
the  original  fiber  and  treated  as  a  one-dimensional  elastic  continuum  was  again  employed.  In  this 
case,  however,  lateral  displacement  of  the  stiffener  was  taken  to  equal  lateral  displacement  in  the 
extended  matrix  along  the  centroidal  axis  of  the  original  fiber  and  Bemoulli-Euler  bending  beam 
theory  was  used  to  describe  the  behavior  of  the  stiffener.  Body-force  field  distributions 
corresponding  to  laterally-loaded  rigid  disks  embedded  in  the  matrix  along  the  axis  of  the  fiber 
were  adopted  as  the  "bond -forces." 

In  the  method  proposed  here,  the  stress  field  in  a  matrix  in  which  a  loaded  fiber  is 
embedded  is  approximated  by  the  stress  field  in  an  identical  matrix  in  which  the  fiber  has  been 
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replaced  by  unknown  distributions  of  point  forces  and  point  dilatations  along  the  fiber's  centroidal 
axis.  The  fiber  is  considered  separately  as  a  one-dimensional  elastic  continuum  whose  stress  and 
deformation  state  is  related  to  matrix  quantities  along  the  fiber-matrix  interface.  The  unknown 
distributions  are  then  solved  for  by  requiring  that  the  fiber  satisfy  equilibrium  and  constitutive 
relations.  Tnis  necessitates  the  numerical  solution  of  a  pair  of  coupled  integral  equations.  In 
contrast  to  the  model  used  by  Muki  and  Sternberg,  this  approach  treats  the  transfer  of  load  between 
the  fiber  and  matrix  in  a  manner  which  explicitly  includes  tangential  tractions  across  the  interface 
and  therefore  affords  one  more  flexibility  in  examining  systems  where  interface  conditions  are  an 
issue.  Furthermore,  the  fundamental  elastostatic  solutions  in  application  in  this  model  are  that  for  a 
point  force  and  a  point  dilatation.  These  solutions  are  much  less  cumbersome  than  the  disk  of 
uniform  loading  (or  laterally-loaded  rigid  disk)  required  in  Muki  and  Sternberg's  model.  These 
factors  make  the  method  presented  below  attractive  for  modeling  a  variety  of  fiber-matrix  systems 
of  interest  to  those  studying  fiber-reinforced  materials. 

1.  LOAD-TRANSFER  TO  AN  ELASTIC  MEDIUM  FROM  AN  INFINITE 

AXIALLY  LOADED  FIBER 

Perhaps  the  best  way  to  present  this  model  is  to  demonstrate  its  application  with  a  simple 
problem,  in  this  case  the  three-dimensional  analog  of  Melan’s  first  problem  from  two-dimensional 
elasticity.  An  infinite  cylindrical  fiber,  with  a  circular  cross-section  of  radius  a,  is  ideally  bonded 
along  its  entire  length  to  a  fully-infinite  matrix  and  subjected  to  a  concentrated  load  F  (see  Fig.  1). 
The  model  is  used  to  solve  for  the  resultant  axial  load  carried  by  the  fiber.  A  cylindrical  coordinate 
system  is  defined  as  shown  in  Fig.  1  with  the  z-axis  coincident  with  the  centroidal  axis  of  the  fiber 
and  the  applied  load  at  the  origin  in  the  negative  z-direction.  Both  the  fiber  and  the  matrix  are 
homogeneous  and  isotropic,  linear  elastic  solids  with  Young’s  modulus  and  Poisson's  ratio  taken 
respectively  to  be  £/and  vy- for  the  fiber  and  Em  and  \m  for  the  matrix. 
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Consider  first  the  fiber  of  the  problem.  In  this  model  the  fiber  is  approximated  as  an 
axisymmetric  elastic  rod  with  a  uniform  axial  stress  c.  This  means  it  is  assumed  that  eq  =  er,  e2, 
Oq  =  cr,  and  c,  are  functions  of  z  only  and  shear  strains  are  ignored.  Under  the  rod  theory 

approximation,  constitutive  relations  for  the  fiber  reduce  to 

a  =  Efez  +  2vfar  (1) 

and 

£/ee  +v/£/ez -(1-2v/)(1  +  v/)ar  =0  .  (2) 

The  fiber,  taken  as  a  free  body,  is  subject  to  a  concentrated  load  F  at  z  =  0  and  to  bonding 
tractions  acting  at  r  -a  between  the  fiber  and  the  matrix.  These  bonding  tractions,  along  with  their 
equivalent  matrix  stresses,  are  a  distributed  shear  stress,  x  =  xj”2(tf,z),  and  a  self-equilibrating 

"pressure",  or  =  o™(a,z).  Throughout  the  remainder  of  this  paper,  field  quantities  in  the  matrix 
will  be  denoted  by  superscribed  ms.  The  rod  is  in  equilibrium  if,  for  all  z. 


7  >  F 

Jifl  o  +  27taJ  X  dz'  =  — sgn(z) 


(3) 


where  the  signum  function  sgn(z)  =  z/|z|.  This  equation  might  of  course  appear  in  different  forms 
depending  on  the  lower  limit  used  in  the  integral  of  the  shear  stress  distribution.  As  written  here, 
however,  (3)  reflects  the  natural  symmetry  of  the  problem.  In  what  follows,  (2)  and  (3)  are  taken 
to  be  the  governing  equations.  Utilizing  (1),  the  governing  equations  contain  three  fiber  quantities 
which  need  to  be  related  to  the  approximate  elastic  field  in  the  matrix  described  below,  namely  eq. 


e2,  and  <jr. 

As  already  stated,  the  elastic  field  in  the  matrix  acted  upon  by  the  loaded  fiber  is 
approximated  by  the  reaction  in  a  fully-infinite  elastic  space  due  to  a  concentrated  force  F  acting  at 
the  origin  in  the  negative  z-direction  along  with  distributions  along  the  z-axis  of  point  forces  and 
point  dilatations,  p(z)  and  q(z)  respectively,  where  the  point  force  distribution  must  be  self- 
equilibrating.  Papkovich  stress  functions,  y  and  <j>,  will  be  used  to  express  this  approximate  elastic 
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field.  In  a  cylindrical  coordinate  system  with  rotational  symmetry  the  expressions  for  radial  and 
axial  displacements  are 


«(r,z)  =  -^^-[zv(r,z)  +  <)(r,z)] 


w(r,z)  =  ^p  z^-v(r)z)-(3-4v)\|r(r,z)  +  ^<J>(r,z)J  .  (5) 

Thus,  using  the  Kelvin  solution  for  a  concentrated  load  in  an  infinite  elastic  space  along  with  that 
for  a  point  dilatation  [8],  the  approximate  elastic  field  can  be  expressed  as 


,z)  +  r~$(r,z)  . 

dZ 


Y(r,z)  =  — ^ - -  J  -  -  -  [F5(0  +  P(0K 

8tt(1  -  vm)  j^p(r,z-QL 


where  p(r,z)  =  Vr2  +  z2  and  8(£)  is  the  Kronecker  delta.  In  this  expression  the  point  forces  are  of 
strength  p(z)  and  act  in  the  negative  z-direction  while  the  point  dilatations  in  expansion  are  of 


magnitude 


-q{i).  Using  these  stress  functions,  expressions  for  all  field  variables  in  the 


matrix  can  be  derived  in  the  form  of  infinite  integrals  of  the  unknown  distributions  multiplied  by 
some  known  difference  kernel,  e.g.  in  the  form 

er<r,z)=  J{A,(r,z-;)[F5(;)  +  pK)]+A3(r,z-C)«K)}^  (7) 


where  the  kernel  functions  Ai  and  A2  are  real  analytic  functions  of  z  for  r  >  0  and  so, 
consequently,  is  e™.  For  those  matrix  quantities  which  shall  be  used  below  in  terms  of  the 
Papkovich  stress  functions  see  the  Appendix.  For  reasons  discussed  in  the  literature  [1]  it  is 
impossible  to  model  the  fiber  as  an  elastic  line  (with  a  =  0)  in  this  problem. 
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The  fiber  quantities  u  and  ar  in  equations  (2)  and  (3)  are  set  equal  to  the  corresponding 
matrix  quantities  at  r  =  a  so  that 

e0=-wm(a,z)  (8) 

a 

and 

ar=o?(a,z)  (9) 

where  we  have  shown  that  each  of  these  quantities  is  analytic  and  must,  therefore,  be  continuous. 
On  the  other  hand,  on  physical  grounds,  and  from  (3),  the  axial  fiber  stress  o  must  certainly  be 
discontinuous.  Then,  assuming  that  Gr  is  continuous,  it  follows  from  (1)  that  ez  is  discontinuous. 
The  difficulty  is  that  the  seemingly  most  natural  expression  for  e2,  that  given  by 

ez=e!?(a,z)  (10) 

in  the  form  (7),  can  not  be  discontinuous  and  so  is  suitable  for  use  only  in  equation  (2).  Some 
alternative  expression  for  ez  in  terms  of  matrix  field  quantities  must  be  adopted  for  equation  (3). 

As  an  indication  that  some  alternative  expression  for  ez  involving  a  discontinuity  would  not 
be  unreasonable,  consider  the  following.  Let  £z  be  the  average,  over  a  disk  of  radius  a  centered  on 
the  z-axis,  of  the  strains  in  the  matrix  ezm(r,z  -  £)  where  quantities  with  a  superscript  asterisk  are 

due  to  a  unit  concentrated  axial  load  at  (r,z)  =  (0,Q.  Calculation  leads  to  an  expression  of  the  form 
£*  =Csgn(z-Q  +  regular  terms  .  (11) 

Except  in  the  range  |z  -  =  O (a)  this  expression  agrees  rather  closely  with  e*m(a,z  -  Q  for  the 

same  load.  However,  this  expression  is  not  suitable  for  the  following  reason.  From  (1)  it  follows 

that  the  "jump"  Aez  in  strain  be  related  to  the  jump  in  az  at  z  =  C,  by  Aez  =  —  Aaz  =  — l —  (for 

Ej  na  Ef 

the  supposed  unit  load).  The  constant  C  in  (11)  does  not  meet  this  requirement.  For  lack  of 
compelling  reasons  for  a  particular  choice  the  effective  axial  strain  in  the  matrix  for  a  unit  axial  load 
at  (0,0  was  taken  to  be 
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♦ 

ez 


(12) 


=  — [sgn(z  -  O' ~  erf  f — 11  +  Cfc. z  ~  0  • 

2na  Efl  \  a  )_ 

A  somewhat  different  expression  will  be  used  in  the  half-space  problem  to  follow.  There  is 

certainly  nothing  unique  about  the  expression  (12)  which  gives  finally  an  expression  for  axial 

strain  to  be  used  in  equations  (1)  and  (3).  Final  numerical  results,  which  agree  very  well  with  those 

in  the  literature,  seem  to  indicate  that  such  results  are  quite  insensitive  to  the  precise  definition  of 
* 

ez- 

Using  the  expressions  established  above  to  relate  the  unknown  fiber  quantities  to  the  elastic 
field  given  by  (6),  and  non-dimensionalizing  with  a  and  F,  the  governing  equations  (2)  and  (3)  can 
be  rewritten  as  a  pair  of  coupled  integral  equations  in  terms  of  the  unknown  distributions  p(z)  and 
<7(z); 

2  fp(?W+  J[r,  1  (r  -  + r,  2  (z  -  =  -r„w  03) 

0  — oo 

oo 

J[r2i(r-;)/Kt)+r22(*-SM?)K  =  -r2,W  04> 

—  oo 

where  the  kernel  functions  rap  (a,P  =  1,2)  are  real  analytic  functions  (see  the  Appendix).  Note 
that  this  system  could  be  solved  analytically  using  Fourier  transform  methods  and  the  convolution 
theorem.  However,  of  concern  here  is  the  establishment  of  methodology  for  more  complicated 
load-diffusion  problems. 

The  system  (13)  and  (14)  is  first  reduced  to  a  set  of  discrete  linear  equations.  Using  the 
symmetry  of  the  distributions,  p(z)  even  and  q(z)  odd,  the  infinite  integrals  can  be  rewritten  in 
semi-infinite  form,  though  without  difference  kernels.  Truncating  infinite  limits  at  appropriately 
large  values  and  approximating  the  integrals  with  a  trapezoidal  quadrature  scheme,  the  two 
equations  are  then  enforced  at  the  discrete  quadrature  points  in  accordance  with  the  Nystrom 
method  [5].  Though  not  truly  a  system  of  first  kind  integral  equations,  (13)  and  (14)  unfortunately 
retain  some  of  the  ill-posed  behavior  inherent  in  all  such  equations.  This  is  dealt  with  by  using 
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singular  value  decomposition  to  solve  the  set  of  linear  equations,  filtering  out  small  length  scale 
instabilities  with  some  unavoidable  degradation  of  the  results  for  small  z  [6].  Under  the 
assumptions  of  this  model,  especially  that  the  fiber  behaves  as  an  elastic  rod  and  has  axial  strain 
given  by  (12),  one  would  not  expect  high  accuracy  near  the  applied  load  in  any  event 

Axial  load  in  the  fiber  is  determined  by  equation  (1)  [or  equivalently  equation  (3)]. 
Rewriting  (1)  in  non-dimensional  form  with  the  now  known  distributions  p(z )  and  q{z)  gives 


^  =  a,w  +  J[a,(z  -  CMC) + a2(,  -  cwopc 


(15) 


where  o0  =  - — j,  Ai  is  discontinuous  at  z  =  0,  and  A2  is  a  real  analytic  function  (see  the 

Appendix).  The  bonding  tractions  across  the  fiber-matrix  interface,  x  and  ar,  can  similarly  be 
determined.  Comparison  between  the  results  obtained  with  this  model  and  those  from  an  exact 
elastostatic  solution  developed  by  Muki  and  Sternberg  [1]  along  with  the  results  of  their 
approximate  model  are  shown  in  Fig,  2-4.  The  results  for  Cf  very  closely  approximate  »he  exact 
solution  and  can  be  shown  to  have  the  same  asymptotic  form  in  the  highest  order  term  as  |z|  -» °°, 
i.e. 


o(z)  ,  ,  Ef 


1- 


(1  +  vm)£//£m+(l  +  V;)(l-2v/) 


sgn(z) 


+  o(z~2)  .  (16) 


The  ratio  of  Young's  moduli  between  the  fiber  and  matrix  is  seen  to  be  much  more  of  an  influential 
factor  than  either  of  the  Poisson's  ratios.  A  comparison  of  cr  in  Fig.  3  points  to  a  shortcoming  of 
the  model  which  gives  a  result  that  is  incorrectly  continuous,  though  it  approaches  the  exact 
solution  asymptotically  for  large  z.  Fig.  4  shows  that  while  the  exact  solution  predicts  a  logarithmic 
singularity  in  t  at  z  =  0,  the  everywhere  bounded  approximate  solution  is  quite  accurate  elsewhere. 
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2.  LOAD-TRANSFER  TO  AN  ELASTIC  HALF-SPACE  FROM  A  SEMI-INFINITE 

AXIALLY  LOADED  FIBER 


Typically  of  greater  interest  in  the  study  of  fiber-reinforced  materials  are  problems  of  load- 
diffusion  from  fibers  in  semi-infinite  half-spaces.  The  problem  of  this  type  solved  here  is  the  three- 
dimensional  analog  of  Reissner's  problem  from  two-dimensional  elasticity.  A  semi-infinite 
cylindrical  fiber,  with  a  circular  cross-section  of  radius  a,  is  ideally  bonded  to  a  semi-infinite 
matrix.  The  fiber  is  normal  to  the  free-surface  of  the  matrix  and  is  subjected  to  a  concentrated  load 
F  away  from  the  matrix  (see  Fig.  5).  A  cylindrical  coordinate  system  is  defined  as  shown  in  Fig.  5 
with  the  z-axis  coincident  with  the  centroidal  axis  of  the  fiber  and  the  matrix  occupying  the  space 
z  >  0.  Both  the  fiber  and  the  matrix  are  homogeneous  and  isotropic,  linear  elastic  solids  with 
Young's  modulus  and  Poisson's  ratio  taken  respectively  to  be  Ef  and  vy  for  the  fiber  and  Em  and 
\m  for  the  matrix. 

In  applying  the  model  to  this  problem  the  procedure  established  above  is  repeated. 
Constitutive  relations  for  the  fiber  under  the  rod  theory  approximation  are  still  given  by  (1)  and  (2) 
and  the  bonding  tractions  along  the  fiber ,  with  their  equivalent  matrix  stresses,  are  denoted  in  the 
same  way.  The  rod  is  in  equilibrium  if,  for  all  z  >  0, 

Z 

na2o  +  2naj x  dz'  -  F  .  (17) 

o 

In  what  follows,  (2)  and  (17)  are  taken  to  be  the  governing  equations  and  £q,  e2,  and  <Jr  are  once 
again  the  fiber  quantities  which  need  to  be  related  to  the  approximate  elastic  field  in  the  matrix 
described  below. 

The  approximate  elastic  field  due  to  a  concentrated  force  F  acting  at  the  origin  in  the 
negative  z-direction  along  with  distributions  along  the  positive  z-axis  of  point  forces  and  point 
dilatations,  p(z)  and  q(z)  respectively,  can  be  expressed  using  the  Mindlin  solution  for  a 
concentrated  load  in  a  semi-infinhe  elastic  half-space  along  with  that  for  a  point  dilatation  [8],  as 
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. |-s«?o!l!sS#*,ii---»i-i,-i“*i'*i*'<''”oi 


(18) 


where  it  is  recalled  that  p (r,z)  =  yjr2  +  z2  .  The  distributions  have  the  same  magnitude  as  before. 
Using  these  stress  functions  results  in  matrix  field  variables  of  the  form 


£?(r,z)  =  J{B1(r,2,0[F5(i;)  +  p(;)]  +  B2(r,z,C)?(Q}^  .  (19) 

0 

While  the  kernels  Bi  and  B2  are  not  difference  kernels  they  are  still  real  analytic  functions  of  z  and 
C,  for  r  >  0  and  so,  consequently,  is  e™ . 

The  fiber-matrix  relations  for  ee  and  Gr  that  were  established  in  (8)  and  (9)  are  still  valid  as 
is  that  for  axial  strain,  ez,  that  was  established  in  (10)  to  be  used  in  equation  (2).  Care  is  used  in 
choosing  an  expression  for  ez  to  be  used  in  (17).  Consider  a  unit  concentrated  axial  load  at  (0,Q.  It 

follows  from  fiber  equilibrium  and  (1)  that  there  must  be  a  jump  in  axial  strain,  Ae*  =  — \ — , 

n  a  Ej- 

across  the  load.  Recall  that  a  superscribed  asterisk  denotes  a  quantity  due  to  a  unit  concentrated 
axial  load.  The  expression  given  in  (12)  satisfies  this  condition,  but  following  the  asymmetric 
nature  of  the  problem  the  effective  axial  strain  in  the  matrix  was  instead  taken  to  be 
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(20) 


where  H(z)  is  the  Heaviside  step  function.  Comparison  with  other  results  in  the  literature  show 
(20)  to  be  an  acceptable  choice. 

Using  the  expressions  established  above  to  relate  the  unknown  fiber  quantities  to  the  elastic 
field  given  by  (18),  and  non-dimensionalizing  with  a  and  F,  the  governing  equations  (2)  and  (17) 
can  be  rewritten  as  a  pair  of  coupled  integral  equations  in  terms  of  the  unknown  distributions  p{z) 


and  q(z); 


]  pW + J[nu(z,Qp(Q+ n,2(z,Q4(0K  =  -n„(r,o) 


0  0 


(21) 


J[n2i(z,C)p(0+nJ2(z.o?(gp;  =  -n21(r,o)  (22) 

0 

where  the  kernel  functions  nap  (a, (3  =  1,2)  are  real  analytic  functions  (see  the  Appendix).  These 
equations  are  solved  for  discrete  values  of  the  distributions  in  the  same  way  as  before  except  that 
the  limits  of  integration  are  already  semi-infinite.  Note  that,  as  previously  alluded  to,  (21)  and  (22) 
do  not  have  difference  kernels  and  could  not,  therefore,  be  solved  using  transform  methods. 

Axial  load  in  the  fiber  is  determined  by  equation  (1)  [or  equivalently  equation  (17)]. 
Rewriting  (1)  in  non-dimensional  form  with  the  now  known  distributions  p(z)  and  q{z)  gives 

^  =  I1(z,0)+J[l1(z,;)p(g  +  Z2(z,CMQ]dC  (23) 

°°  0 

p 

where  a0  = — j,  Zi(z,g  is  discontinuous  at  :  =  and  I2  is  a  real  analytic  function  (see  the 
na 

Appendix).  Results  from  (23)  are  compared  with  those  from  Muki  and  Sternberg’s  model  [2]  and 
shown  in  Fig.  6. 
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CONCLUDING  REMARKS 


Comparison  with  an  exact  solution  for  the  three-dimensional  analog  of  Melan's  first 
problem  [1]  shows  that  the  fiber  load-diffusion  model  gives  good  results  for  both  axial  load  in  the 
fiber  and  tangential  tractions  on  the  fiber-matrix  interface  at  distances  from  the  applied  load  greater 
than  approximately  one  fiber  radius.  The  discontinuous  nature  of  the  normal  interface  tractions, 
however,  is  not  adequately  accounted  for.  To  demonstrate  the  model's  application  to  fiber  load- 
diffusion  problems  involving  half-spaces  the  three-dimensional  analog  to  Reissner’s  problem  was 
examined.  The  results  were  compared  to  those  from  Muki  and  Sternberg's  approximation  [2,3] 
and  found  to  be  in  agreement.  It  is  hoped  that  the  model  will  prove  useful  in  the  study  of  other 
more  complex  fiber-matrix  systems  which  more  closely  resemble  those  observed  in  actual  fiber- 
reinforced  materials.  In  particular  it  is  anticipated  that  the  model  will  be  useful  in  studying  systems 
with  prescribed  interface  conditions  other  than  ideal  bonding. 
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APPENDIX 


Listed  below  are  those  matrix  quantities  that  are  needed  for  the  implementation  of  the  model 
in  this  paper.  They  are  given  in  terms  of  the  Papkovich  stress  functions  which  in  turn  are  given,  in 
terms  of  the  unknown  distributions  p(z)  and  q(z),  by  equation  (6)  for  the  full-space  problem  and 
equation  (18)  for  the  half-space  problem. 

L>,z)  =  i^|-[ry(r,2)  +  <|>(r,z)]  (Al) 

r  hm  or 


1  +  v 


-  2(1  -  2vm)— \|/(r,z)  +  ^-j(|)(r,z) 


dz 


dz 


(A2) 


<  (r,  z)  =  ^j[z\|/(r,  z)  +  <t >(r,  z)]  -  2vm  Jj z) 


(A3) 


1 G (r, z') dz'  =  |- [z>(r, z')  +  <f)(r, z'fl  *  ■ -  2(1  -  Vm ) ] ^{r,z')dz’ 


d_ 

dr 


(A4) 


0  0 
Note  that,  as  stated  in  the  text,  each  of  these  quantities  is  a  real  analytic  function  of  z  for  r  >  0. 


The  three-dimensional  analog  of  the  Melon  problem. 

The  kernel  functions  for  the  pair  of  coupled  integral  equations  (13)  and  (14)  are  arrived  at 
by  substituting  (Al)  -  (A4)  [with  \j/  and  0  given  by  (6)]  into  the  governing  equations  (2)  and  (3) 
using  the  fiber-matrix  relations  established  in  the  body  of  the  paper.  To  non-dimensionalize  with  a 

and  F  the  following  substitutions  are  made:  z->az,  pa-»ap,  8(Q-»-5(£), 

a 

p(Q  ->  --p{Q,  and  q(Q  Fq(Q.  Letting  pz  =  Vl  +  z2  ,  the  kernel  functions  are 
a 

r,i(z)  =  -erf(z)-» - 1-  (Xj  — j  +  0-2  ~ t  (A5) 

Pz  Pz  Pz 
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r12  (z)  -  a3  — 3  +  a4  — y 

Pz  Pz 


r2i(z)  -  a5  —  +  a6  — 
Pz  Pz 

„  ,  v_  1  1 

^22  (z)  -  a7  — 5" +  a8  ~ J 
Pz  Pz 


where 


«i 


=  (l  +  v  )ii£-l+-v/(1  _}Vm) 


Em  2(1 -Vm) 


an  = 


2(l-vm)L  '  2 


v 


~(l  +  vm)-X 


a3  =  l  +  v/-(l  +  vm)-^- 


«4  =-(l  +  vm)-^-3v/ 


«5=0  +  vm) 


2(1 -vm) 


-2v 


f 


Ef  (l-2v/)(l  +  v/)(l-2vm) 


2(1 -v.) 


«6  = 


_ 


(1  + V  1 

1  T  vm 

E,  )  3(l-2Vy)(l  +  V/) 

U-O 

Em  2(1  -v„) 

a7=(l-2v/)  1  +  Vy -(l  +  vm)— ■ - 


a8  =  -3v/(l  +  vm)-^-3(l-2v/)(l  +  v/)  . 


The  kernel  functions  in  the  expression  for  axial  load  in  the  fiber  (15)  are 


Ai(z)  =  sgn(z)  -  erf(z)  +  Pi  -y  +  P2  "T 

Pz  Pz 


(A6) 

(A7) 

(A8) 

(A9) 

(A10) 

(All) 

(A  12) 

(A  13) 

(A  14) 

(A15) 

(A  16) 

(A  17) 
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(A  18) 


A2(z)  ~  P3T  +  P4- T 

Pz  Pz 


where 


g  |  v  \Ef  v/(1  2v™) 

P‘~(  V")*-  2(1  -v.) 


(A  19) 


p2  =  ai  =  — - -  V/- — (l  +  v_)— 

2  *  2(l-vm)[/  2^  m)Em 


(A20) 


p3  =  V/-(l  +  vm)-£ 


(A21) 


3  Ef 

P4=a4=-(l  +  vm)-f-3v/ 


(A22) 


The  three-dimensional  analog  of  Reissner's  problem. 

The  kernel  functions  for  the  pair  of  coupled  integral  equations  (21)  and  (22)  are  arrived  at 
by  substituting  (Al)  -  (A4)  [with  \jr  and  <))  given  by  (18)]  into  the  governing  equations  (2)  and  (17) 
using  the  fiber-matrix  relations  established  in  the  body  of  the  paper.  The  same  substitutions  are 
made  to  non-dimensionalize  with  a  and  F.  The  kernel  functions  are 


n„(*.C)  =  -H(i-c)erf(2-c)+ifl±^')|^o1(2.c)+ 


+  4(l-vm)“3(2’C) 


(A23) 


n!2(z’0  =  7O  +  vm)7Tfl4(z.0+  "^'^5(z-0  +  T^6(z<0 


(A  24) 
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.  _  v  f  ( 1  +  v  ^E,  ,  x  (l  +  Vf)(l-2v,-) 

n21W)  =  f  ^  /n,(r,Q-L_  /A  "q2(z.C) 

8  V1_VWfcm  8(1- VmJ 


1 


(A25) 


n22(z-0  -  -j-(!  +  vm)-^Q4(z,Q  -  U\  +  Vyjfl  -  2v/)q5(z,Q 


+  T(1  +  vm)7T-f28(z.C) 


(A26) 


where  Q,(z,Q  (i  =  1...8)  are  continuous,  analytic  functions  given  below.  Letting  p  =  ^1  +  (z  -£)2 
andpsJl  +  (z  +  02, 


r.,(..Q  ~30^</  +  ^  3(l-4v„-4K)(z  +  5)+6z 


4(l-4vm  +  2v;L)(2  +  0+4z  3(2-5)  4(l-vm)(2-C) 


(A27) 


/./■)_  30zC(z  +  Q  6z5(2  +  5)+12v„z-9(z-0 
fl2(2'a - ?  p3 - 


,65-(l-2v„)(3-4vm)(z  +  0 . 3(z-5)  (l-2vJ(z-5) 
+  +_^ - - 


+4(1-v„)(1-2v„) 


(z  +  «3+5(2  +  Op 


(z  +  C+pfp3 


„  -6z?(z  +  C)  (3-4vm)z  +  C  z-C  2(1 -v„)  2(1 -v„) 

^z3vz>^/  —5  —3  3  /  ,  r  .  — /  *•  , 


(A28) 


P  (r  +  ;  +  p)p  (z-;  +  p)p 


+4(1 -vm) 


(A  29) 


17 


-30z(z  +  g  ,  12z(z  +  g-3+ 12vm  .  2 - 8vm  ,  3  2 


cum)*  t;?  '+ 


"—5 - * — r - t 

P  P  P 


(A  30) 


^  /  r\...30z(2  +  0  6z(z  +  C)  +  9  3  +  4vm  3  1 

q5(z*?)— T - 7 — +~7~~7  7 


O6(*.0* 


_  -6z(z  +  0  1  1 


"^V 


0  (  r\  _  -64(z  +  0  (3-4vm)(z-C)  z-C  4(l-vm)(l-2vm) 

- r? - =5 - 3~  +  — ;  -T-x- — 


(z  +  C  +  p)p 


(A31) 


(A32) 


(A33) 


flg(*.C)  = 


_  -6z(z  +  Q  3  -  4vm  1 


+  7T  • 


(A  34) 


The  kernel  functions  in  the  expression  for  axial  load  in  the  fiber  (23)  are 


Z1(2,g  =  H(2-g[l-erf(2-y]  +  l(|^)|^n1(z,0  +  --^p;n2(z,0  (A35) 

®  V1  ~  Vm  /  ‘-‘m  4l1-Vm) 


S:(-.0 = 7(1  +  + ^-n5(2,5) 


(A36) 


where  i2,-(z,£)  (i  =  1,2, 4, 5)  are  as  defined  above. 
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Fig.  1  -  Infinite  cylindrical  fiber  embedded  in  an  infinite  elastic  matrix. 


0.0  1.0  2.0  3.0  4.0  5.0 

z/a 


Fig.  2  -  Load-diffusion  results  in  the  three-dimensional  analog  of  Melan's  first  problem  for  the 
model  presented  in  this  paper  and  Muki  and  Sternberg's  exact  and  approximate  formulations. 
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Fig.  3  -  Normal  stress  across  the  fiber-matrix  interface  in  the  three-dimensional  analog  of  Melan’s 
first  problem  for  the  model  presented  in  this  paper  and  the  exact  solution. 


z/a 

Fig.  4  -  Shear  stress  across  the  fiber-matrix  interface  in  the  three-dimensional  analog  of  Melan's 
first  problem  for  the  model  presented  in  this  paper  and  the  exact  solution. 
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